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Abstract. - We study a one-dimensional extended Hubbard model with longer-range Coulomb 
interactions at quarter-fiUing in the strong coupling limit. We find two different charge-ordered 
ground states as the strength of the longer range interactions is varied. At lower energies, 
these CO states drive two different spin-ordered ground states. A variety of response functions 
computed here bear a remarkable resemblance to recent experimental observations for organic 
TMTSF systems, and so we propose that these systems are proximate to a QCP associated 
with T = Q charge order. For a ladder system relevant to Sr\4,Cu2i04,\, we find in-chain CO, 
rung-dimer, and orbital antiferromagnetic ordered phases with varying interchain couplings and 
superconductivity with hole-doping. 



Electron crystallization, or charge ordering (CO) due to interactions, is an issue of enduring 
interest in condensed matter physics. The study of the conditions favoring CO, along with 
its competition with metallic and/or superconducting states constitutes a problem of wide- 
ranging interest for a host of real systems [1] composed of weakly coupled single chains /ladders. 

In this work, we study this issue within an extended quarter-filled Hubbard model on a 
linear chain, described by, 

i^cff ~tY^{ClCi+i,„ + h.c) + {U~ 2zP)Y^ n.^mi +VY^ n^n.+i + P ^1 ^^'"^+2 

i,fT i i i 

In one-dimension, the spin fluctuations are those of an ideal S = 1/2 XXX AF chain, while 
the charge fluctuations are described by the Hamiltonian 

He = -t ^(c|c,+i + h.c) + {V - J 14) n,n,+i + n,n,+2 (2) 

i i i 

that describes a model with frustrating interactions. In ID, the projected fermions are spinless 
fermions with a hard- core constraint. This model has been considered as a model for studying 

(*) E-mail: M. S .Laadaiboro . ac .uk 
(**) E-mail: sanjaylOthp.uni-koeln.de 

© EDP Sciences 



2 



EUROPHYSICS LETTERS 



the effect of frustration on electron crystallization [2]. For narrow-band systems, we consider 
the limit t << (F — J/4), P. In this regime, we employ an extension of the trick used for the Id 
next-nearest neighbor Ising chain: for t = 0,we notice that with {V — J/4) > 2P, the ground 
state is the usual CDW (Wigner) crystal for n = 1. With 2P > {V — J/ 4,), however, the 
dimerized state (Peierls) is the ground state, written schematically as (11001100....). Splitting 
this in a slightly different way, we have [...(01)(10)(01)(10)...]. Associating a pseudospin 
T = 1/2 operator, with = +1 for (10) and -1 for (01), the state is antiferromagnetic and 
doubly degenerate in terms of the rf. For small t, this is an attractive trick because (in spin 
language) the transverse term does flip the r?, but cannot break a pair. So one obtains, 

ifeff = - E[2*^r + - J/^ - '^pytn^i] (3) 

I 

This is just the Ising model in a transverse field, which has been studied extensively in ID. 
If (y — J/4— 2P) < 0, the ground state is ferromagnetically ordered in r^, i.e, it corresponds to 
a Wigner CDW. For (y— J/4 — 2P) > 0, the Pcicrls dimcr order results in the groimd state. At 
(V — J/4 — 2P) < 2t, the quantum disordered phase has short-ranged pseudospin correlations, 
and is a charge "valence-bond" liquid. The quantum critical point at {V — J/4 — 2P) = 2t 
separating these phases is a dcconfincd phase with gaplcss pseudospin (r) excitations, and 
power-law fall-off in the pseudospin-pseudospin correlation functions. Correspondingly, the 
density-density correlation function has a power-law singular behavior at low energy, with an 
exponent a = 1/4 characteristic of the 2D Ising model at criticality. For P = 0, the metallic 
phase for 1/ < 2f is a Luttinger liquid, and in this limit, the low-energy physics is qualitatively 
similar to that of the usual t — J model. The "Mott" insulating state for V > 2t has Wigner 
CO in the groimd state, and the M-I transition is of the Kostcrlitz-Thoulcss type [3]. 

The full Hamiltonian in our case for the strong-coupling limit is now given by 

ifeff =-I^[2irf + {V~-2P)tM+i] + Si-Si+i (4) 

I I 

To study the magnetic phases, we adapt the Ogata-Shiba [4] technique for our case. This 
is possible if J << t,V, in which case, the pseudospin part is first solved exactly (this is 
possible because of the known exact solution of the ID transverse field Ising model), and 
the exchange part is then treated as a perturbation. Writing the total wavefunction as a 
product of a spin and pseudospin wavefunction (where the spin wavefunction is defined in a 
Hilbert space of dimension 2^), i.e, 1-0 >= |t > ®\S >, and following standard degenerate 
perturbation theory, the spin degeneracy is lifted by the correction (of order 1/L): 

< H,e >'= -2t<T^>+Y, Jm+i(Si.Si+i - 1/4) (5) 

; 

where the average < .. >' denotes that the average is taken over the exact ground state |r > of 
the pseudospin part above, i.e. < A >'=< t\A\t > and Jij+i = "^(^ ^ -f ^ 2P) < rfr^^^^ >. 

An interesting fact now emerges: Wigner CO (FM order of t) results in an HAFM S = 
1/2 spin model with the Hamiltonian Hg = J^^S,.S,+i. This gives rise to a gaplcss AF 
gromid state for the spin degrees of freedom. The charge (holon) excitations arc gapped: this 
corresponds to a linear confining potential for holons. On the other hand, Peierls dimerization 
in the charge sector (AF Neel order of r) gives rise to dimerization in the spin sector, with 
the Hamiltonian = JEJl + (-l)*(5]Sj.Sj+i. 

Translated into fermion variables, this yields a sine-Gordon problem with = 2n, and 
describes an instability to a singlet pinned ground state commensurate with the Peierls CO 
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setting in at higher energies. The elementary excitations are soUtons carrying = ±1. 
ScaUng theory predicts a dimer gap, Ad ~ 6"^^^. Exactly at = 2tt, the SG model has just 
two 5^ = breather excitations with opposite parity [5], the lowest, even parity breather 
being degenerate with the /S^ = ±1 soliton doublet, forming a S* = 1 triplet, while the second 
odd-parity breather is a singlet with a gap, ^/SA^- It is important to notice that both charge 
and spin order arise from long range Coulomb interactions, and do not involve an electron 
phonon coupling mechanism. 

Let us consider the implications of having the CO state in the high-T regime, where one 
could imagine the system to be effectively one-dimensional. In particular, we want to look 
at the LJ,T dependence of the various response functions at high-T. Using the exact solution 
of the pseudospin model in ID, the high T (in the "quantum critical" regime) behavior can 
be explicitly derived [6]. In fact, near Ising criticality, the response function, x{^) — r~^l^ 
where r = (x^ -|-t^)^/^ (with the velocity v set to unity). This relation is still valid away from 
criticality in the "short range" region, r « A^-, where is the pseudospin (in our case 
charge gap) gap of the ID-TFIM. Using this asymptotic form, we have 



Xcrit(0,a;) = -^^^B\A -iS,l- 2A) (6) 



where S = -^j^, and A = 1/16 is the conformal dimension. B{x,y) is the beta function. 
In the quantum critical region, an illuminating form is 

. x(o,o) 

1 — ZW/rR, + K^^^ — (w/Wi)^ 

where Tr = (2tan(7r/16)fcBT//i)e-'^-/'==^'^, wi = 0.795(A:BT/fi) and ^ =l.28{ch/kBT)e^-/''^'^ , 
are determined solely by T and the fundamental natural constants, as expected in the QC 
regime. Here, A^- is the energy gap to charge excitations in the Wigner/Peierls CO states de- 
scribed above. This represents the collective charge susceptibility, and the optical conductivity 
follows directly from (7{u)) = — ia;x(0, w), giving. 



x(o,o) ^ 

Tr (1 - a;Vw?)2 + (w/Fr) 



(^H = — ..2/. ,2X2 . /../-n ^2 («) 



The corresponding frequency-dependent dielectric function is obtained from e(a;) = 1 -|- 
{'inia{ui)/ui), and the electronic contribution to the Raman scattering is estimated therefrom 
to be given by Ib.{u}) = lm(l/e(0,a;)), for light polarized along the chain axis. In terms of the 
charge susceptibility, this is simply, 

^ e{Lo) i + (l£^m)2^2(^^2^) ^ ^ 

where F(^,T) = ^,_^y^.y.+^^/r^y. ■ 

X" {k, oj) has its maximum value at Wm = wi — i{ujf/Tn), implying that the collective mode 
broadens and shifts to higher energy linearly in T with increasing T at high temperatures. 
Further, the T-dependent damping rate of the collective mode correlates well with the relax- 
ational peak seen in transport, underlying their common origin. In fact, the dc resistivity is 
linear in T at high T, with "insulating" features showing up at lower T . In our picture, these 
are collective (longitudinal) bosonic charge-density modes in the high-T quantum critical re- 
gion above an incipient CO transition (expected to occur at low T). In fig.(l), we show the 
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Fig. 1 - A three-dimensional plot of the Raman intensity /^j versus a scaled energy uj/T and temper- 
ature T for parameter values of the original model which constitute a gap = COSfes. 

electronic Raman lineshape as a function of lo/T. The sharp low energy peak corresponds to 
the collective charge density fluctuation mode of the CO ground state. 

The corresponding frequency-dependent dielectric constant also shows an explicit w/T 
scaling in the QC regime, or generally, at high-T, it shows strong T-dependence. From fig. (2), 
we see that it becomes w-independent at high T, but appreciably increases as T is lowered, 
with a maximum at a; ~ T. 




Fig. 2 - A three-dimensional plot of a(Lo)/uj versus a scaled energy uj/T and temperature T for 
At = O.OSfca. The limiting behaviors are: (1) <j{uj) ~ ^"^y^ for high T, and ~ for low T. 
(2) e{Lo,T) = 1 + ^ for high T, and e'{Lo,T) ~ ^ for low T. 

The fact that organic charge transfer salts [7] exhibit features very similar to those found 
above has interesting implications. In light of our results, these anomalous features can now 
be identified with proximity to an underlying quantum critical point associated with charge 
(Wigner/Peierls) ordering. We recall that very recent work [8,9] shows that the dimerized 
insulating state in TMTSF systems has charge order at low T. Interestingly, e'(0,a;) indeed 
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shows appreciable increase as T is lowered, further supporting an interpretation based on 
proximity to an underlying CO ground state [7]. Hence, we conclude that observation of 
these features in TMTSF systems constitutes strong evidence that this system is close to a 
putative QCP associated with charge order. Observation of dimerized/Neel ordered AFM 
states co-existing CO states at low T are also naturally understood in light of the analysis 
above [9]. 

Two-chain Ladders We now consider the strong coupling version of a coupled two-chain 
ladder system, with each chain being described by H as in eq.(l). In the strong coupling 
limit, where each chain is described by a TFIM for charge degrees of freedom, the coupled 
chain model is constructed as follows. For U ^ oo, and V,P>t (but {V — 2P) comparable 
to t), the charge degrees of freedom of the fermionic problem for each chain are described by 
an effective pseudospin model on n-n bonds, via the effective Hamiltonian, 

jjchain ^ _ ^[2fr/ + {V- J/4 - 2P) t/t/+i] (10) 

3 

Rotating the pseudospin axis such that t'^ ^ , ^ —t^ and coupling two such chains via 
an interaction coupling U±_ and a two-electron interchain transfer t_i_, we have the effective 
Hamiltonian for the charge sector of the two chain system as 

H = -Y,[^tTl, + {V-J|^-2P)Tl,T^+^J- J2 {U±rlarl,+ t±{rl,Tl,+Tl^Tl,)], {11) 

j,a j,a,bjia 

where a,b = 1,2 is the chain index. Denote the in-chain pseudospin coupling as J = {V ~ 
J/4 — 2P) and the inter-chain pseudospin coupling as J± = U±. Here, we study the strong 
coupling version of this problem in two limits (see below). The weak-coupling problem is 
studied elsewhere [10]. 

For the case of | Jx| >> | J|, t± , the 2 chain system can be better thought of as strongly- 
coupled rungs which are weakly coupled to their neighboring rungs. Thus, we treat J as 
a weak perturbation on the zeroth-order system of rungs defined by the large coupling Jj_, 
giving Hef f = Ho + Hi where 

Ho = -hY,Tl, + J^ ^laTlb 
j.a j,a,b^a 

Hi = -jJ2^!,arUi,a-YT.(4,ar-, + ii.c) (12) 

j,a j,a,bjta 

where the effective magnetic field is given hy h = 2t > 0. 

LEH for J_L < . For J_L < and h « J±, we find that the triplet state |-|-) = -^(1 ti 

) + I it)) a-nd the singlet state |— ) = -^(| t|) — | arc degenerate on any rung and are 

separated from all other states by a large gap of order J±. Thus, these two states define the 
subspace which will determine the low-energy physics of the system. Identifying a pseudospin- 
1/2 operator with the low-energy subspace on each rung, we treat the Hamiltonian Hi as 
a perturbation (to second order in J/J±) to obtain the LEH as 

H - 2^1- ^ij2 4/^2 ^3+1 -^2^^3 

3 3 
j2 l2 t2 j2 _ 27)2 
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We find that t± acts as the strength of a Zeeman-sphtting hke term in the LEH. Bosonising 
this, we obtain a sine-Gordon Hamiltonian with a cosine potential in the dual (9) field and a 
magnetic-field term 

H = lm<l>f + id.ef] - £^ cosp^e - ^d,<i> (14) 

We note that bosonisation of the general XYZ Hamiltonian results in the appearance of 
an additional 4fcir Umklapp term [11], cos/32<?i', which is irrelevant for a finite t± and is 
hence ignored in what follows. When t± is below a certain critical value, incommensurate 
Wigner charge order (ordering of the £,^) occurs [12]. Above this critical value, a spin-flop 
transition orders the system in the x direction (i.e., ordering of the ^'^) via a Kosterlitz- 
Thouless transition. For < Stt, the cosine in the dual field is a relevant perturbation and 
orders the dual field. The magnetic-field term ex t± leads to a ground state with charges which 
are coherently delocalised on the diagonals of each pair of nearest-neighbor rungs; this is an 
orbital antiferromagnet-type ordering with circulating currents in plaquettcs [3, 11]. 

LEH for J_L > . For Jj_ > 0, and h > 0, the triplet state \+) = | tt) is the low 
energy state on any rung. For h = 0, we find that the triplet states |-|-) (defined above) and 
|— ) = I II) are degenerate. Thus, we can again identify these two states as the subspacc which 
determines the low-energy physics of the system. For h << J, we again identify a pseudospin- 
1/2 operator with the low-energy subspace on each rung, and treat the Hamiltonian Hi as 
a perturbation (to second order in J/J±) to obtain the LEH as 

^ = -ijE^Kj+i-2iE^|- (15) 
j j 

This is just the ID TFIM (with ferromagnetic Ising coupling). In the ordered phase, the 
ground state has in-chain Wigner CO and dimers on every alternate rung. The disordered 
phase is a gapped, short-ranged charge-dimer liquid. At < = J^/4Jx, the quantum critical 
point describes a gapless charge-dimer liquid with Lu/T,vk/T QC scaling, exactly as was 
described before. Transposing the results obtained before, we conclude that the dc resistivity, 
optical conductivity, electronic Raman and dielectric responses will be exactly described by 
the same scaling functions (eqs.(6)-(9)) with the gap, A^, now being the CO gap of the ladder 
problem {H in eq.(15)). Very interestingly, exactly such behavior is observed in undoped 
ladder system 5ri4Cu2404i [13] and attributed to a longitudinal, collective charge fiuctuation 
mode, exactly as described here. 

LEH for hole-doped ladder. Upon doping the ladder with holes, while a single hole ex- 
periences a linear confining potential in the Wigner (Ising-like) or Peierls (dimerized) CO 
bac;kground, a pair of holes on the same rung is free to propagate. One can then describe 
the hole-pair as a hard-core boson, representing its creation and annihilation operators using 
the spin- 1/2 operators cr^; the local charge density is then described by a^. Following [3], we 
find the LEH describing the dynamics of such hole-pairs to be the XXZ model in an external 
magnetic field 

j 

where th ^ J'^ / J± is the pair- hopping matrix clement, u\i is the Coulomb interaction between 
pairs on nearest-neighbour rungs and /x is the chemical potential of the holes. The phase 
diagram of this model is known [3]: for = and Uh > th, the ground state is an insulating 
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CDW of hole pairs. Beyond a critical /Xc = f{uh, th), the system has a ground state described 
by Bosc condensation of hole pairs. In fact, from the bosonisation analysis of the equivalent 
5=1/2 XXZ model in an external Zeeman field, we know that < crfaf^j. >~ r~^/" and 
< o't'^i+r where a = 1/2 — tt"-^ sin~^(2uh/ih). Clearly, for a < 1, the ground state 

has dominant superconducting correlations. This is true for both the cases described above: 
in the first case, wc have a Bose condensate of intrachain pairs of holes, while in the second 
hole pairs on individual rungs Bose condense, describing two possible superconducting types 
in the ladder system. This finding matches our conclusions obtained from a weak coupling 
analysis [10], and thus constitutes a generic feature of undoped/doped strongly correlated 
ladder systems. 

To conclude, we have explored the strong-coupling limit of strongly correlated single chain 
and two-leg ladder models using a variety of methods. Our results strongly suggest that 
prototype examples like organics (TMTSF) and Srij^Cu2AOii lie in close proximity to an un- 
derlying QCP associated with charge order, and constitutes an advance in our understanding 
of their physical responses in a new theoretical framework. 

* * * 
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